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Abstract 

We consider the operator H := idt + V ■ (cV) in an unbounded strip Q in R 2 , where c(x,y) £ C 3 (f2). 
We prove an adapted global Carleman estimate and an energy estimate for this operator. Using these 
estimates, we give a stability result for the diffusion coefficient c(x,y). 

AMS 2000 subject classification: 35J10, 35R30. 

1 Introduction 

Let (1 = R x (-|, |) be an unbounded strip of R 2 with a fixed width d. We will consider the Schrodinger 
equation 

C Hq:=id t q + V ■ (c(x,y)Vq) = in Q = fix (0,T), 
(1.1) I q(x,y,t) =b(x,y,t) on E = 3fi X (0,T), 

I q(x,y,0) = qo(x,y) on fi, 

where c(x,y) € C 3 (fi) and c(x,y) > c m i„ > 0. Moreover, we assume that c and all its derivatives up 
to order three are bounded. If we assume that go belongs to H 4 (Q.) and 6 is sufficiently regular (e.g. 
b e fl' 1 (O,T,.H'8 +e (0fi)) n H 2 (0,T,H% +e (dfl)) and some additional conditions), then dTTTJ admits a 
solution in f/ 1 (0, T, Hi +£ (Sl)). We will use this regularity result later. The aim of this paper is to 
give a stability and uniqueness result for the coefficient c(x, y) using global Carleman estimates and 
energy estimates. We denote by v the outward unit normal to fi on T — <9fi. We denote V = T + U 
where F + = {(x,y) 6 T; y = | } and F~ = {(x,y) £ T; y = — |}. We use the following notations 
V ■ (cVm) = d x (cd x u) + dy(cdyu), Vu • Vt> = d x ud x v + d y ud y v, d u u = Vu ■ ^. 

We shall use the following notations Q = Q x (0, T), Q = SI x (-T, T), E = T x (0, T), E = T x (-T, T), 
A(-Ri) := {$ 6 L°°(n),0 < Ri < \\®\\ L ac (n) }, and A(R 2 ) := {$ £ L°°(n), ||*|U<«»(o) < ^2}, where R ± 
and R2 are positive constants with R\ < i?2- 

Our problem can be stated as follows: 

Is it possible to determine the coefficient c(x,y) from the measurement of d I/ (dtq) on T + ? 

Let q (resp. q) be a solution of (| 1 . 1 f) associated with (c, 6, qo) (resp. (c, 6, qo)) satisfying some regularity 
properties: 

• dtq, V{dtq) and A(9 t q) are in A(i?2), 
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• qa is a real valued function in C 3 (f2), 

• go and all its derivatives up to order three are in A(i?2) • 
Our main result is 



|c - c\ 2 H i {n) < C\dv(d t q) - d v (dtq) 



i ! ((o,T)xr+)i 



where C is a positive constant which depends on (Q, F, T, Ri,R2) and where the above norms are weighted 
Sobolev norms. 

The major novelty of this paper is to give an H l stability estimate for the diffusion coefficient with only 
one observation in an unbounded domain. 

We prove an adapted global Carleman estimate and an energy estimate for the operator H with a 
boundary term on T + . Such energy estimate has been proved in [23] for the Schrodinger operator in a 
bounded domain in order to obtain a controllability result. Then using these estimates and following 
the method developed by Imanuvilov, Isakov and Yamamoto for the Lame system in [16], [17], we 
give a stability and uniqueness result for the diffusion coefficient c(x,y). Note that this stability result 
corresponds to a stability result for three linked coefficients (c, d y c and d y c) with only one observation. 
For independent coefficients, in our knowledge, there is no stability result with one observation. 
The method of Carleman estimates was introduced in the field of inverse problems in the works of 
Bukhgeim and Klibanov (see [T], [3J, [IS], [20]). The first stability result for a multidimensional inverse 
problem (for a hyperbolic equation) was obtained by Puel and Yamamoto [24] using a modification of 
the idea of [3J. 

For the non stationnary Schrodinger equation, [2] gives a stability result for the potential in a bounded 
domain. For the stationnary Schrodinger equation, we can cite recent results concerning uniqueness for 
the potential from partial Cauchy data (see for exemple [TH] and the references herein). 
In unbounded domains Carleman estimate with an internal observation has been proved for the heat 
equation in [6]. 

A physical background could be the characterization of the diffusion coefficient for a strip in geophysics. 
Indeed if we look for time harmonic solutions of (|l.ip . the problem can be written, after some changes 
of variables as the reconstruction of a non local potential P in a strip for the operator —A + P. Few 
results for inverse problems exist in a two-dimensional strip (see [§]). For the layer R n x [0, h] with n > 2, 
several results exist for the stationnary inverse problems (see [5], [10], [8], [13], [15], [25], ...). 
On the other hand, we can link our problem to the determination of the curvature function for a curved 
quantum guide (see [T2], [TJ, [TT], ...). 

This paper is organized as follows. In section 2, we give an adapted global Carleman estimate for the 
operator H . In section 3, we prove an energy estimate and we give a stability result for the diffusion 
coefficient c. 

2 Global Carleman Estimate 

Let c = c(x,y) be a bounded positive function in C 3 (ft) such that 

Assumption 2.1. c(x, y) £ A(iZi), c and all its derivatives up to order three are in A(i?2). 

Let q — q(x, y, t) be a function equals to zero on dQ, x (— T, T) and solution of the Schrodinger equation 



We prove here a global Carleman-type estimate for q with_a single observation acting on a part F + of the 
boundary T in the right-hand side of the estimate. Let (3 be a C 4 (Q) positive function such that there 
exists positive constant C pc which satisfies 

Assumption 2.2. • |Vj9| € A(Ri), <9„/3 < on T~ , 

• /3 and all its derivatives up to order four are in A(i?2). 



idtq + V ■(c(x,y)Vq) = f. 



• 2K(D 2 /3(C, 0) - cVc • V/3|C| 2 + 2c 2 |V/3 • C| 2 > C pc |C! 2 , for all C G C 



where 
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Note that the last assertion of Assumption 12.21 expresses the pseudo-convexity condition for the function 
j3. This Assumption imposes restrictive conditions for the choice of the functions c(x, y) in connection 
with the function /3. Note that there exists functions satisfying such Assumptions; indeed, if we consider 



c(x,y) £ < / £ C (fi); 3ro positive constant, 



-/ d y fdyl3 > r > 0, 

/ dyfdyf3((^) 2 + 1) + 2f(dy y f3 + {OyPf) > T > 0. 

then a function /3(x,y) = /3(y) is available (for example, c(x,y) — {j^t + l)e~ y and /3(x,y) = e y ). 

Similar restrictive conditions have been highlighted for the hyperbolic case in [5T], [25] and for the 
Schrodinger operator in [14] 

Then, we define /3 = /? + K with K = m||/3||oo and m > 1. For A > and t £ (— T, T), we define the 
following weight functions 

(2-1) = „(x,y,t) = {T + t)(T _ t y 

Let H be the operator defined by 

(2.2) Hq:= id t q + V ■ (c(x,y)Vq) in Q = fix(-T,T). 

We set tp = e~ sv q, Mtfi = e~ sr, H(e sv ip) for s > and we introduce the following operators 

(2.3) MiV := idtip + V ■ (cVV>) + s 2 c\Vr 1 \ 2 tp, 

(2.4) M 2 i/> := isd t r)ip + 2csS7n ■ Vip + sV ■ (cVr/)?/>. 
Then the following result holds. 

Theorem 2.3. £et _ff, Mi, M 2 be the operators defined respectively by \2.2\) . (2.3$ . \2.1$ . We assume 
that Assumvtions \2. il and \2.2\ are satisfied. Then there exist Ao > 0, so > and a positive constant 
C — C(Q, r, T, Cpc, Ri, R2) such that, for any A > Ao and any s > so, the next inequality holds: 

(2.5) s 3 A 4 j T j e- 2s,1 <p A \q\ 2 dx dy dt + s\ j T j e" 2s >| Vq\ 2 dx dy dt + ||Mi (e~ s "q) 

r-T r pT 



+ \\M 2 (e- s \)\\ 2 L2(d) <C 

-t Jr+ 



sX / / e 2sv ifild^qf d v (3 da dt + I e 2sr> \Hq\ 2 dx dy dt 
J —T Jt+ J —t J n 



for all q satisfying Hq £ L 2 (Q x (— T, T)), q £ L 2 (-T, T; ffj(fl)), <9„ 9 £ L 2 (—T, T; L 2 {Y)). 
Proof. If we set ?/> = e~ sv q, we calculate = e~ 3V H(e sv '^) and we obtain: 

Mip = + M 2 ?/> 
with Mi and M 2 defined respectively by (|2.3|l and (|2.4|) . Then 



(2.6) f [ \Mip\ 2 dxdydt = f f \Mup\ 2 dx dy dt + [ [ \M 2 ip\ 2 dx dy dt 

J-TJn J-TJn J-rJn 

+ 23i( f [ Mnp M^ dx dy dt), 
J -t Jn 



where z is the conjugate of 2, K (z) its real part and (z) its imaginary part. We have to compute the 
scalar product in (|2.6[) 



r-T 



5R( / / Mnp M 2 ip dx dy dt) = I 1V + I 12 , + + I 2V + I 22 , + I 2Z , + I 31 , + I i2 , + I : 



XV 
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Then, we have 

(■T 



(2.7) I IV = ®-(yJ J {idtip)(-isd t r) ip) dx dy dtj = -| J J dttv \ip\ 2 dx dy dt. 

I l2 , = K (2is J J c dtipVr) ■ Vi/> dx dy dt^j 

= sQ (^J J c ' ^ t% P dx dy dt^j — sQ (^J J c Vr; • Vip dt4> dx dy dt 

After an integration by parts with respect to the space variable in the first integral and to the time 
variable in the second integral, we obtain 

(2.8X12' = -s3 (J J V • (cVt?) V dt^ dx dy d?J + (J J cip V{d t rj) ■ dx dy d?J . 

(2.9) I 13 > = ^(f J V • (cVr?) ip dti> dx dy dtj . 
Note that I 13 / vanishes with the first term of 7i2' . 

(2.10) hi' = U (-is J J <9t7? ^ V ■ {cVip)dx dy dt^j 



— sxs 



T 

cip S/(d t rj) ■ V-0 dx dy dt 



t Jn 



Using integrations by parts, we obtain 

(2.11) I 22 > = 2sJ? (^j j cV ■ (cVV>)V?7 ■ V^)dx dy dt\ 



= -s\ [ f cp\Vip\ 2 (y ■ (c 2 V/3) + \c\VP\ 2 ) dx dy dt 
J -t Jn 

+ sf f cd v T) |9„V| 2 da dt + 2sA 2 f f ipc 2 \S7 f3 ■ Vip\ 2 dx dy dt 

+ 2s ak j i i P c Yl °- ' ; ' dx d y dt 

= -s\ f [ <p\ViP\ 2 (V ■ (c 2 V/3) + \c 2 \Vf3\ 2 ) dx dy dt 

J -T Jn 

- sA f [ c 2 tp d u (3 \d„ip\ 2 da dt + 2s\ 2 f f ipc 2 \X7f3 ■ Vip\ 2 dx dy dt 
J -t J an J-tJci 

+ 2sAK ( f [ (pD 2 f3(X7ip, V^) dx dy dt 



4 



Using integration by parts, we obtain 

(2.12) 7 23 / = sSR j V ■ (cW) V • (cVt?) ^ da; dy di 

= sX [ f ap \ViP\ 2 (V ■ (c\7f3) + \c\Vf3\ 2 ) dx dy dt 



t Jn 

/ <^V ■ (cV ■ (cV/3)V/?))H 2 d:r dy dt 

f y>c|V/3| 2 V ■ (cV/3)\4>\ 2 dx dy dt 
t Jn 

4 rT r „\3 /-T 



/ / iy3c 2 |V/3| 4 |V)| 2 da; dy dt - — f f tpV ■ (c 2 |V/3| 2 V/3)|-0| 2 dx dy dt 
J-rJn 2 J- T Jn 



TJB ' J-TJSl 

f-T r 

/ y>cV/3 ■ V(V ■ (cV/3) + Ac|V/3| 2 )|-0| 2 ds dy di 

T Jf2 

/ ipV ■ (cV(V ■ (cV/3) + Ac|V/3| 2 ))|V>| 2 dz dy dt. 
-t Jn 



sX 
2 

sA 
~2 

sA 
~2 

sA 
~2 

sA 
T 

And we obviously have 

(2.13) I 31 ,=s 3 $t([ I c(-id tV nV V \ z ipdxdydt) =0. 



(2.14) J 32 < = 2s 3 5t^y J c 2 \Vr)\ 2 ipVr)-Vip dx dy dt 

= s 3 / / c 2 |Vt7| 2 V77 • V|V| 2 dx dy dt 

= s 3 A 3 f f op 3 (|V/3| 2 V -{cVP) + VI3-V{c\Vf3\ 2 ))\4)\ 2 dx dy dt 
J -t Jn 

+ 3s 3 A 4 f f c 2 ip 3 \X7f3\ 4 \ip\ 2 dx dy dt. 
J -t Jn 

(2.15) I 33 / = s 3 K (J j f c |V77|VV ■ (cVr?)^ dx dy dtj 

= -s 3 A 3 f f c <^ 3 |V/3| 2 V • (cS7/3)\iP\ 2 dx dy dt 
J -t Jn 



cV|V/3| 4 |V| 2 ds; dy dt. 



So, by (p71) - (j2~T5)) . we get 

rT r \ fT 



K (/ / M2 ^ d:r dj/ dt ) = 2sA2 / / VPc2 ' v,s ' v ^\ 2dx d y dt 

- sA f f <p cVc-V/3\Vip\ 2 dx dy dt 

J -t Jn 

+ 2sX$t (J J (pD 2 f3(Vifj, Vi/j) dx dy d?J 

+ 2s 3 A 4 f [ c 2 ip 3 \Vf3\ 4 \ifj\ 2 dx dy dt 
J —t Jn 

— sX / c 2 tpdvfl\dv^\ 2 da dt + X, 

J -t J an 



where 



2 

sX 
2 

sA 3 
2 



X = — - / / dttV M 2 dx dy dt + 2s9 / / c^V{d t r)) ■ V^dx dy dt 



t Jn 



ipV ■ (cV ■ (cV/3)V/3))|^| 2 di dy dt 
<pc\VI3\ 2 V ■ (cV(3)\4>\ 2 dx dy dt 



IT f I V c2 \VP\ 4 M 2 dxdydt-^- f 



■1± 

2 

*A 
T 



t Jn J —T J fi 

r T ,■ 

\ <pcX7/3 ■ V(V • (cV/3) + Ac|V/3| 2 )H 2 d:r dy dt 

T 

ipV ■ (cV(V ■ (cV/3) + Ac|V/3| 2 ))|V>| 2 dx dy dt 



T J n 



+ s 3 X 3 [ I ^ 3 cV/3- V(c|V / 3| 2 )|V| 2 dz dy dt. 



Recall that: 



\\MiI>\\ lH q) = ||M 1 V|| 2 2( Q ) + ||M 2 ^|| 2 2((5) + 2K(Af 1 ^M 2 ^). 



Then : 

(2.16) ||A/n/>|| 2 2(c5) + ||M 2 VC 2C q) + / T / ^ c2 |V/3 ■ W>| 2 dz dy dt 



t Jn 



rT r \ rT 



+4sA» ^ y <pD 2 (3{S7ip, V4))dx dy dtj ~ 2sX J J V c Vc ■ V/3|V-0| 2 ds dy dt 
+4s 3 A 4 y J c 2 V > 3 \Vf3\ 4 \i>\ 2 dx dy dt-2s\ [ [ c 2 <pd„l3\d„Tp\ 2 da dt <\\MtP\\ 2 l2i5] +2\X\ 



ll£2(Q) 

-rJn J -t J an 



Taking into account 

• |/3| + |V/3| + | V(V • (cV/8))| + |V • (V(V • (cV/3)))| < C(n, T, T, i? 2 ) in SI, 

• |ftti7| < C(TV, < C(T)^ 2 , p < C(2> 3 , ^ 2 < C(7V, 

• \s5{[ f cTpV{d t rj) -V^dx dy dt)\ <C{T)s\ f f c<p\V f3 ■ Vip\ 2 dx dy dt 

J-TJn J-TJn 

+C(T)s\ f f cip 3 \ip\ 2 dx dy dt, 
J-t Jn 

where C(Sl, T, T, R2) is a positive constant depending upon Q,T,T, R2 and C(T) is a positive constant 
depending upon T. Therefore we obtain the following estimation for X: 



\X\ < C(Sl, r, T, R2) 



(sX 4 + s 3 X 3 )f [ <p 3 \tP\ 2 dx dy dt + s\ [ f tp\V(3 -Vij}\ 2 dx dy dt 
J-TJn J-TJn 



The two terms of the previous estimate of \X\ are neglectable with respect to 

s 3 A 4 f [ c 2 Lp 3 \Vf3\ 4 \Tp\ 2 dx dy dt or sA 2 f f ipc 2 \V/3 ■ Vip\ 2 dx dy dt, 
J-TJn J-TJn 



for s and A sufficiently large. Using Assumption I2.2I we have 



4sAK ^ y <pD 2 f3(Vip, Vip)dx dy dtj - 2sA j j <p c Vc • S//3\S/^\ 2 dx dy dt 
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+4sA 2 f [ ipc 2 \Vf3 ■ Vipfdx dy dt > C pc s\ [ [ cp\Vip\ 2 dx dy dt, 
so (|2.16p becomes 

l|AMll» ( Q) + \\ M ^\\l*(Q) +4s 3 \ 4 J T ^J^c 2 ^\X7(3\*\i}\ 2 dx dy dt 

rT 

|2 j„ j„. jj. o„\ / / 2 a „/.|2 



+sA / (p\Vip\ dx dy dt < 2s\ / c 1 ipd v fi\d v i>\ da dt + \\Mip\\l 2(d) . 
J-tJu J-TJdn 

Recall that <9„,3 < on c(x,y) £ \(Ri) n A(J? 2 ), |V/3| £ A(Ri) and V = e~ s,| <7, then the proof is 
complete. 

□ 

3 Inverse Problem 

In this section, we establish a stability inequality and deduce a uniqueness result for the coefficient c. 
The Carleman estimate (|2.5[1 proved in section 2 will be the key ingredient in the proof of such a stability 
estimate. 

Let q be solution of 

( idtq + V ■ (cVg) = in Ox(0,T), 

(3.1) ^ q(x,y,t) = b(x,y,t) on 9fix(0,T), 

and 5 be solution of 

r i<9 t q + V ■ (cV§) = in Qx(0,T), 

(3.2) i q(x,y,t) = b(x,y,t) on 9fix(0,T), 

I q(x,y,0) = go(z,2/) in 



where c and c both satisfy Assumption 12.11 If we set u = q — q, v = ftu and 7 = c — c, then u and v 
satisfy respectively 

( id t u + V ■ (cVu) = V ■ (7V5) in fix(0,T), 

(3.3) { u(as,y,t) = on aflx(0,T), 
[ y, 0) = in fi, 

f id t v + V ■ (cVu) = V ■ (-yVdtq) = f in !)x(0,T), 

(3.4) ^ u(a;,y,t)=0 on 9flx(0,T), 

Assumption 3.1. go is a real valued function in C 3 (Q) 

We extend the function v on Q x (— T, T) by the formula v(a;, y, t) = — «(a;, y, —t) for every (x, y, t) £ 
f2 x (— T, 0). Note that this extension is available if the initial data is a real valued function. For a pure 
imaginary initial data, the right extension is v(x,y,t) = v(x,y,—t). Note that these extensions satisfy 
the previous Carleman estimate. 

3.1 Energy Estimate 

We assume throughout this section that 7 6 Hq(Q). We introduce 

(3.5) E(t)= / e- 23v(x - v ^\dtu(x,y,t)\ 2 dx dy + [ tp-\x, y, t) (,-*•">&<*>*) \d t Vu{x, y, t) \ 2 dx dy. 

Jn Jn 

In this section, we will give an estimation of E(0). 
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First Step: We first give an estimation of J n e- 2sv( - x ' yfi) \dtu(x,y,0)\ 2 dx dy. 
We set ip = e~ sv v. With the operator 

(3.6) Mn/> = id t ^ + V-(cVV0 + s 2 |V7?| 2 ^, 

we introduce, following [2J, 

J=2q(J J Miiplp dx dy dt 

Assumption 3.2. dtq, V(dtq), A(dtq) are in A(B,2). 

We have the following estimate 

Lemma 3.3. We assume that Assumvtion \3 ,S\ is satisfied. Then there exists a positive constant C 
C(fi, r, T, Ri, R2) such that for any A > Ao and s > so, we have 

1 = f e- 2svix - y - o) \dtu(x,y y 0)\ 2 dx dy 
and 

\X\ < Cs~ 3/2 \- 2 f e- 23 ^ y ' 0) (\j\ 2 + \V~f\ 2 ) dx dy 

+Cs- 1/2 \- 1 f T f e~ 2ar] ip d v /3 \d„v\ 2 da dt. 

J-T JT+ 

Proof. In a first step, we calculate T 

X = 2 ^(y J {idtip^ + V ■ (cV^)? + s 2 |Vr7| 2 V V>) dx dy dt\ 

= 23? (^J J dt iM> dx dy dl 

° f ,2 

/ dt\4>\ dx dy dt 

-t Jn 



\ip(x,y,0)\ 2 dx dy 
Jn 



So, we have 



X= f e- 2s,li ^ yfi) \d t u(x,y,0)\ 2 dx dy. 
In a second step, we estimate I. Using Young inequality we can write 

\I\ < J \M^\ 2 dx dy dt^j 2 (^J J e~ 2sv \v\ 2 dx dy di 

< C(T)s~i\- 2 (\\Mii> f L3 + s 3 A 4 J J e~ 2sn ip z \vf dx dy dt 

with C(T) a positive constant which depends on T. Then with the Carleman estimate (|2,5|l proved in 
section 2 we have 

f e~ 2sr > 'ipd„f3\d n vf da dt+ [ f e~ 2s7) |V ■ (jVd t q) dx dy dt] , 
-T J r+ J-tJsi ) 

where C = C(fl, F + , T, Ri, R2) is a positive constant. Using Assumption 13.21 since 

e -2B V {x,y,t) < e -2«,(a,»,0) for a „ t e (_ T;T ) j 



we obtain for s and A sufficiently large the estimate 

(3.7) \X\ < C s- 3/2 \- 2 f e - 2s "(^^>°)(| 7 | 2 + |V 7 | 2 ) dx dy 



in 

+C s- 112 ^ 1 [ I e" 2s V d v /3 \d u vf da dt. 

J — T JT+ 

where C = C(Q,, T + , T, Ri , R2) is a positive constant. □ 

Second Step: We then give an estimate of J n ip^ 1 (x, y, 0) e- 2s ^ x ' v ' 0) \d t S7u(x,y,0)\ 2 dx dy. 
We denote 

(3.8) E(t):= [ c^-\x,y,t) e- 2 ^^\Vv{x,y,t)\ 2 dx dy, 



where tp 1 = — . We give an estimate for E(0) in Theorem 13.51 In a first step we prove the following 
lemma : 



Lemma 3.4. Let v be solution of \3.4\ l in the following class 

v £ Ci^T],!! 1 ^)), d v v 6 L 2 {Q,T,L 2 {T)). 
Then the following identity holds true 

e~ 2srl f ip~ 1 dtv dx dy dt 



-2s 

forfsHlipL). 



E(t) -E(k) = -25ft ( / / 

\Jk Jfl 

-Sft (J J c e _2s "(-4sA + 2\ip~ 1 )d t v V/3 • V« dx dy dt 

f c e-^ip-^t^X/vl 2 dx dy dt + f f ce- 2av d t (v~ 1 )\Vv\' 2 
Jn Jn Jn 



Proof. Since v is solution of (3.4) note that d t v — if — iV • (cVw). Therefore, we obtain the two following 
equalities. 

(3.9) U [ f e~ 2sv f tp^dtV dx dy dt = $t ( -i [ [ e' 2av f l P ~ 1 V ■ (c VtJ) dx dy dt) , 

Jk Jn V J k, Jn J 

(3.10) Sft (J I ce~ 2srl dtv (-4sA + 2A ¥ ?~ 1 )V / 3- dx dy dt) 

= 5ft (i J J c e' 2sv J (-4sA + 2\tp~ 1 )S7(3 ■ Vv dx dy dt\ 

-Sft (i f f c e _2s "(-4sA + 2Avp _1 )V ■ (c\7v)X7f3 ■ Vv dx dy dt 
\ J K Jn 



We multiply the first equation of (|3.4p by e 2sr> tp 1 d t v and we integrate on (k, t) x Q., where [k, t] C [0, T\ . 
So, if we consider the real part of the obtained equality, we have 

= Cft j j e~ 2s,) tp'^dtvf dx dy dt) = 5ft (J j e~ 2s " f tp^dtv dx dy dtj 

-5ft (J J e~ 2er > ^ _1 V ■ (cVv)d t v dx dy dt 
Then by integration by parts, we obtain 

= Eftf^ J e~ 2a,) f tp^div dx dy dtj j c\7v ■X7(e~ 2sri ip^dtv) d. 



= » / / e- 2sv f ip^dtv dx dy dt) - 2sK / / c ^V 28 " ^uVu • Vr? dx dy dt 
\J K Jn ) \J K Jn 

-AS^y f c e~ 2sv ip^dtvVv -V/3 dx dy dt^j + ^(^f f c e~ 2sv tp'^v V (d t v) dx dy dt 

Note that 

E{t) - E{k) = f c d t {e~ 2sv v^lVvl 2 ) dx dy dt. 

Therefore we have 

= 5ft ( / / e~ 2s,) / <p~ 1 d t v dx dy dt) - 2s5R (/ [ c tp' 1 e' 2sv d t vVv • Vn dx dy dt 
\J K Jn J \J K Jn 

-A5R ( f [ c e~ 2sv V ' 1 d t vVv ■ V/3 dx dy dt) + -E(r) - -E(k) 



<> 



2 v ' 2 



I c e 2sv tp 1 d t n\\7v\ 2 dx dy dt — i f ce 2sv d t (ip 1 )|Vu| 2 dx dy dt, 



and the proof of Lemma 13,41 is complete. □ 
Theorem 3.5. Let v be solution of \3.$ in the following class 

v e c([o,T],H\n)), d v v e L 2 (o,r,L 2 (r)). 



We assume that Assumvtions \27l\ and \27w are checked. Then there exists a positive constant 
C = C(tt,r,T,Ri,R 2 ) > such that 



(3.11) E(0) < C 



[ e~ 2s \ d u j3 |cU| 2 do dt + sX [ [ e~ 2s,1 \f\ 2 
■t Jr+ J J q 



for s and X sufficiently large. 

Proof. We apply Lemma [3 . 4 1 with k — and t — T. Since E(T) — 0, we obtain 

E(0) = 2$t(f f e' 2 " 1 f tp^dtV dx dy dt\ ~^(f f c e~ 2s "(-4sA + 2X(p~ 1 )d t v X7f3 ■ Vv dx dy d?J 

+2s / / c e~ 2sr '(p~ 1 dtn\Vv\ 2 dx dy dt — / / ce~ 2s,, <9t(y> _1 )|Vi;| 2 . 
Jo Jn Jo Jn 

We give now estimates of the four integrals in the previous equality. 

First integral: B := 25ft (7 Q T / n e'^ip' 1 / d t v dx dy dt) . 
Using (|3.9|l . we have : 



j-T i- \ / r T 

(3.12) B = 25ft ' 



e 2sr 'ip 1 / dtv dx dy dtj = 5ft (^f f e 2 " v ip 1 f d t v dx dy dt^j 
+5R (-if f ^""ip^f V • (c W) dx dy d?J . 



Recall that if we set V = e~ sr >v, then Mip = e~ sn '//(e 3 "^) = Miip + M 2 ip for s > with 

MiV := idtip + V ■ (cVVO + s 2 c| Vn\ 2 i/j, 
M 2 ip := isdtnip + 2csVrj ■ V?/> + sV ■ (cVri)ip. 

So (|3.12[) becomes 

B = 25ftQ| f e~ aOT V -1 / dtv dx dy dt^j =^(f f e~ 2sv ip' 1 f e sv (sd t nrp + dtii) dx dy di 
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+5R 



Note that 
Then, we obtain 



e" 21 " 7 ^" 1 / e s? '(sV ■ (cVr])4> + c s 2 \Vrj\' 2 ip + 2c s V77 ■ Vip + V ■ (cV^)) da; dy di 



3? 



e s > _ 7 e BV (idtip + V ■ (cVV)) dx dy dt 



e arl ip f e Brl (sdtr] — i sV ■ (cVr?) — i c s \Vrj\ )ip dx dy dt 



-Si 



e 23,7 <p / e s? '2c s V77 • dx dy dt 



idtip + V ■ (cVVO = - s 2 c|Vt?| 2 V' 



B = K 



+5R 



e 2s7 'y> x / e s, (Mit/) — s 2 c| V/7| 2 t/>) da; dy dt 

Jn 

e -2 SV(p -if e s*i( s9tr] _ i sV . ( cVr? ) _ j c s^Vryl 2 )^ da; dy dt 

T 



e 2sv tp 1 f e" v 2c s V77 • Vip dx dy dt 



'0 Jfi 

If we come back to the function v, the previous equality becomes : 



B = R[-i 



-0? 



(/e S? V 1 Mi(e s,| i;) — s z c\Vrj\^e ^ sr, tp 1 / u) da; dy dt 



2 — 2srj -1 



e V 9 / (sdtfj — i sV ■ (eV/7) + i c s |V??| )i> dx dy dt 

-2sti -1 

V 

Then there exists a positive constant C = C(r2, F, T, .R2) such that 



+K ( -i / / e / 2c s V?? ■ Vu da; dy dt 

Jo in 



(3.13) 



e S, V / ft 1 ' dx dy dt\ < C 



-2s?) I pi2 



f\ J dx dy dt + WA'hie-^v)]] 2 ^ 



Q 



+s 3 \ 4 



J J e 2sv <p 2 \v\ 2 dx dy dt + s\ J J e 2arl \SIv\ 2 dx dy dt 



Second integral: D := 5R { / Q / n c (-4sA + 2Av3 _1 ) e^ 2s?) (l + tp- 1 )d t v V/3 • Vv da; dy dt 
We denote by p := — 4sA + 2A( ( 9~ 1 . Using (|3.10[1 . we have 



2D = U 



c e- zsr >pd t v V/3 • Vu dx dy dt) + / / c e SI7 / pV/3 • Vv da; dy dt 

in 



-5R 



If we introduce tp = e sv v, we get: 



c e 2av pV ■ (cV?J)V/3 • V« da; dy dt 



2D = R 



-dt i 



'0 Jn 

■T 

'0 Jn 



ce~ 2B "pVv ■ V/3(sd t rje sv ^ + e sv d t ip) dx dy dt ) +0? ( i / y ce -28 " Vw ■ V/3 / dx dy dt 



ce'^pVv ■ X7(3[s 2 af)\X7r l \ 2 e l3V + 2sce'" 7 V^ • Vr? + se^V • (cVr?) + e s "V • (cV^/j)] 
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Therefore 



2D = 3t(-ij J ce~ ar >pVv-Vf3[M 1 il> + M2'ip] dx dy dt^ J Jce- 2av pV', 

Thus there exists a positive constant C = C(Q, F, T, Ri, R2) such that 

(3.14) \2D\ <C s\ J J [|Mi(e- s "u)| 2 + |M 2 (e- s "w)| 2 ] dx dy dt 

+s\J J e- 2av \Vv\ 2 + sxj J e- 2a "\f\ 2 dx dy dt 

Two last integrals: 

There exists a positive constant C — C(f2, F, T, R2) such that 

(3.15) I f f ce~ 2sv dt(<p~ 1 )\Vv\ 2 dx dy dt\ < C f f e~ 2sv \Vv\ 2 dx dy dt 

Jo Jn J Jq 

and 
(3.16) 



77 

Jo Jn 



cdtrje 2sr, (p 1 \Vv\ 2 dx dy dt\ < C s / ipe 2sr| |Vi>| 2 dx dy dt. 



7/ < 



Using now the Carleman estimate of Theorem 12.31 and Lemma 13.41 from (3.9)-(3.12), we deduce the 
existence of a positive constant C = C(Q,, F, T, Ri, R2) such that: 



e 2 " v ifi d v fi \d v v^\ da dt + sA / / e 2sr> \f\ 2 dx dy dt 
-t Jr+ J Jq 



□ 



E(0) < C 
and the proof is complete 

3.2 Stability Estimate 

Now following an idea developed in [16] for Lame system in bounded domains, we give an underestimate 
for E(0). We adapt the proof of lemma 3.2 of [TB] to an unbounded domain. 

Assumption 3.6. • go and all its derivatives up to order three are in A(i?2) 

• \V/3-Vq \ eA(Ri) 



Lemma 3.7. We consider the first order partial differential operator 

(P g)(x,y) = d x q (x,y)d x g(x,y) + d y (x,y)d y g(x,y),P g := Vq ■ Vg 

where qo satisfies Assumptions lff.il Iff. 61 Then there exist positive constants Ai > 0, si > and C = 
C(Q, F, T, Ri, R2) such that for all A > Ai and s> s\ 

s 2 A 2 / ipo e~ 2sno \g\ 2 dx dy < C I ip^ 1 e~ 2sr]0 -Po5| 2 dx dy 



with rio{x,y) := r/(x,y,0), <po{x,y) := tp(x,y,0) and for g G Hq(SI). 

Proof. Let g G H&(Cl). We denote by w = e" s?,0 p with r) := 7](x,y,0) and Q w = e" s? '° Po{e svo w), so 
we get Qow — swPor/o + Pow. Therefore, we have: 



(p Qow Qow dx dy = s tfo \ w \ l-Fb^o | dx dy + / ip \Pow\ dx dy 
n Jn Jn 



+2sK 



01 ^ 



w Porjo Pqw dx dy 
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= s 2 A 2 / ipo \w\ 2 (Vqo ■ V/3) 2 dx dy + / ip 1 \Pow\ 2 dx dy — sA I Vgo ■ V/3 Vgo • V(|«j| 2 ) da; dy. 
Jn Jn Jn 

So, integrating by parts, we obtain 

/ 'Pa 1 Qow Qow dx dy = s 2 A 2 / ipo \w\ 2 (Vqo ■ V/3) 2 dx dy + / ipQ 1 \Pow\ 2 dx dy 
Jn Jn Jn 

+sA / V ■ (Vg ■ Vf3Vq )\w\ 2 dx dy. 
Jn 

Thus 

/ Vo 1 e~ 2avo \Pog\ 2 dxdy>s 2 \ 2 f ip \Vf3-Vq \ 2 e- 2s,10 \g\ 2 dx dy + sX f V- {P /3 Vg )e- 2s? '° \g\ 2 dx dy. 
Jn Jn Jn 

Using Assumptions 13, ll and 13.61 we can conclude for s and A sufficiently large. □ 

Then, we deduce the following result. 

Lemma 3.8. Let u be solution of \3.4ty - We assume that Assumvtions \2.2\ Iff, 1\ and 1 3. 61 are satisfied. 
Then there exists a positive constant C = C(fi, T, T, Ri, R2) such that for s and A sufficiently large, the 
two following estimates hold true 

(3.17) s 2 A 2 f ip e- 2s,)0 |7| 2 dx dy<C I ipo 1 e- 2am '\d t u{x,y,Q)\ 2 dx dy, 

Jn Jn 

(3.18) s 2 A 2 / e~ 2s?)0 |V 7 | 2 dxdy<C [ ip^ 1 e' 2 " 10 (\V{d t u(x, y, 0))| 2 + | 7 | 2 ) dx dy, 

Jn Jn 

/or 7 eH 2 (fi). 

Proof. We apply Lemma l3.7l to the first order partial differential equations satisfied by 

• 7 given by the initial condition in (|3.4p 

Po7 : = d x qod x 7 + d y qod y "f = id t u(x, y, 0) - 7 Ago, 

• d x y given by the x-derivative of the initial condition in (|3.4[) 

Pod x y ~ d x q d x (d x y) + d y q d y (d x j) 

= id t {d x u(x,y,Q)) - d x "f(Aq + d xx qo) - d y jd xy qo -jd x (Aq ), 

• d y y given by the y-derivative of the initial condition in (|3.4p 

Pod y j := d x q d x (d y "f) + dyqodyidy-y) 

= id t (d y u(x,y,0)) - d y j(Aq + d yy qo) - d x jd xy q - jd y (Aq ). 

Then using Lemma [XT] and Assumptions 13 . ll 13 . 6 1 the proof of Lemma 13.81 is complete. □ 

Theorem 3.9. Let q and q be solutions of L'J. 1\) and H3.2\) such that c — c G H (Q). We assume 
that Assumptions \2. 1\ \2.2i \3.2l \3. 1\ and \3. 61 are satisfied. Then there exists a positive constant C = 
C(fi, r, T, Ri, R2) such that for s and A sufficiently large, 

[ ipo e~ 23V0 (\c-c\ 2 + |V(c-c)| 2 ) dx dy < C f [ <p e- 2sr ' 'd v f3 \d v {d t q - d t q)\ 2 do dt. 
Jn J-tJt+ 
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Proof. Adding (|3.17[) and (|3.18[) we obtain using the estimate (|3.7|l for \T\ and the energy estimate (|3.1ip 
for E(0) 

s 2 X 2 f <p e- 2s ™ (jV 7 | 2 + | 7 | 2 ) dxdy<C f ^ e' 2 ^' (|V(ftu(x, y, 0))| 2 + \8 t u(x, y, 0)| 2 ) dx dy 
Jn Jn 

<C(\X\+E(0)) 

< Cs-' i/2 \- 2 [ e- 23,)0 (| 7 | 2 + |V 7 | 2 ) dx dy + Cs- 1 ' 2 ^ 1 [ [ e _2 *V d v f3 \d v v\ 2 da dt 
Jn J-tJt+ 



+Cs 2 \ 2 



'■ f T [ e- 2s \ d u P \d v v\ 2 do dt + CsX [ [ e- 2sn \f\ 2 . 
J -t Jr+ J J q 

So we get 

s 2 A 2 f ipo e- 2sm (|V 7 | 2 + | 7 | 2 ) dx dy < Cs 2 \ 2 C f e~ 2sv <p d v /3 \d v v\ 2 da dt 



t Jr+ 



+Cs\ J J e"" s " |V ■ ( 7 Va t g)| 2 dx dy dt 

< Cs 2 \ 2 J J e~ 2s " cp 8,13 \d„v\ 2 da dt + Cs\ j j e~ 2s " (|V 7 | 2 + | 7 | 2 ) dx dy dt. 
Then, for s and A sufficiently large, the theorem is proved. □ 

Remark 3.10. This result is also available for the heat equation in bounded or unbounded domains. 
Note that all the previous results proved mix (-|, |) are available in R n x ( — |, |) for n > 2 if we 
adapt the regularity properties of the initial and boundary conditions. 
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